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Abstract 

The new phenomenon of semiquantum chaos is analyzed in a clas- 
sically regular double- well oscillator model. Here it arises from a 
doubling of the number of effectively classical degrees of freedom, 
which are nonlinearly coupled in a Gaussian variational approximation 
(TDHF) to full quantum mechanics. The resulting first-order nondis- 
sipative autonomous flow system shows energy dependent transitions 
between regular behavior and semiquantum chaos, which we monitor 
by Poincare sections and a suitable frequency correlation function re- 
lated to the density matrix. We discuss the general importance of this 
new form of deterministic chaos and point out the necessity to study 
open (dissipative) quantum systems, in order to observe it experimen- 
tally. 
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1 Introduction 



Recently the concept of semiquantum chaos has been introduced in order to 
characterize a particular form of deterministic chaos [|l|. Namely, dynami- 
cal systems with both quantum and classical degrees of freedom may show 
irregular behavior due to their generic coupled nonlinearities 0, ^ ^ . The ap- 
pearance of both quantum and classical degrees of freedom (d.o.f.) generally 
can be thought of as a more or less valid approximation, depending on the 
physical circumstances, describing a truly complex nondissipative quantum 
system. 

For example, the authors of Ref. [0] considered the zero momentum 
(long wavelength) part of the problem of pair production of charged scalar 
particles by a strong external electric field. In this limit the problem could 
be reduced to a classical oscillator interacting with a quantum mechani- 
cal one through a biquadratic coupling; i.e. a two d.o.f. system, one clas- 
sical and the other purely quantum. They represent the electromagnetic 
vector potential and the charged matter field, respectively. In this case 
and for a physically motivated subset of all possible initial conditions, the 
dynamics can be mapped onto an equivalent completely classical problem 
with two effective d.o.f. Their motion is described by two coupled non- 
linear second-order differential equations, which yield regular or irregular 
trajectories depending on initial conditions and suitably chosen model pa- 
rameters. Clearly, an approximation to a purely quantum system, a quan- 
tum field theory in fact, has been treated here and eventually gives rise to 
deterministic chaos in an equivalent low- dimensional classical system. 

The latter observation presents the starting point of our investigation of 
semiquantum chaos in the one- dimensional double-well problem. It can be 
considered as the zero- dimensional limit of the scalar Higgs field model 
To study the entropy production and thermalization in strong interactions at 
ultrarelativistic energies a related scalar quantum field theory has recently 
been treated in the time- dependent Hartree-Fock approximation (TDHF) 
equivalent to a Gaussian variational principle ||^, (see also further refer- 
ences therein). The resulting field equations in the long wavelength limit, 
i.e. neglecting spatial variations, reduce to a coupled pair of Hartree-Fock 
type equations for one-dimensional double-well oscillators. They represent a 
semiclassical approximation to the full quantum problem in terms of spatially 
homogenous fast and slow modes, respectively. The equations for a single 
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such oscillator will be given in Section II and their solutions studied in detail 
in the following. 

At this point it seems sufficient to state that mean-field type TDHF 
quantum corrections to the classical equation for the double-well oscilla- 
tor, which behaves regularly in the classical limit, introduce characteristic 
nonlinear features including additional effectively classical d.o.f. They, how- 
ever, describe quantum fluctuations of the basic coordinates and conjugate 
momenta. Such additional d.o.f. are an essential element to allow for the 
possibility of semiquantum chaos. We recall that a classical one-dimensional 
anharmonic oscillator, which can be described by a set of two autonomous 
flrst-order differential equations, cannot have chaotic trajectories according 
to the Poincare-Bendixson theorem [|I|, which rules out chaotic flow in a 
bounded region in two-dimensional (phase) space. 

The time-dependent variational principle |J7[ using generalized Gaussian 
trial wave functions has, in fact, already been applied in quantum chem- 
istry for quite a while (see Refs. P] for some early and more recent work, 
respectively, citing numerous further relevant references). However, it has 
only recently been pointed out that the parameters (parameter functions in 
field theory 0) specifying the Gaussian wave functions (functionals in field 
theory) are governed by a set of equations according to the variational prin- 
ciple which describe an equivalent classical Hamiltonian system 0. It is this 
equivalent classical system which may show what has been termed semiquan- 
tum chaos, even if the subset of equations belonging to the proper classical 
limit (neglecting quantum fiuctuation d.o.f.) produces only strictly regular 
motion. Furthermore, it is generally believed that the full quantum dynamics 
obeying the Schrodinger wave equation, which is a linear partial differential 
equation (however, a functional l.p.d. equation in field theory), is always 
regular (see Refs. p, |[ for further discussions on this point). Thus, the 
question of quantum chaos usually pertains to more subtle aspects of the dy- 
namics which refiect some traces of deterministic chaos in the corresponding 
classical system 0, |^ , if at all. 

Presently, we want to address the question instead, rather independently 
of the behavior of the full quantum dynamics or even a regular classical limit 
thereof, whether or how the onset of semiquantum chaos is related to the 
breakdown of the underlying semiclassical approximation. Correspondingly, 
one may reconsider the often posed problem of whether or where in the 
relevant classical phase space there are regions in which the system can safely 
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be or even has to be studied using a semiclassical approach |T0 

Finally, we want to draw attention to Ref. [Q. There, the possibility 
of chaotic behavior without classical counterpart is pointed out in a many- 
body system undergoing multiple resonant tunneling. More specifically, an 
electron cloud in a three-well heterostructure was represented by an effective 
one-particle wave function obeying a nonlinear Hartree equation. Clearly, the 
inherently semiclassical description gives rise to the essential nonlinearity of 
the second-order partial differential equation again. However, the authors of 
this work speculate that the many-body character of the system, i.e. infinitely 
many d.o.f. already on the classical field level, might stabilize the reported 
semi quantum chaos in this case to persist even in a fully quantum treatment. 
This necessitates a time-dependent approach to quantum field theory, which 
has not been worked out for the nonrelativistic problem at hand to date. We 
will come back to this in the final section where we discuss one possible first 
step beyond our present simple double-well oscillator model in the direction 
of a quantum field theory, which is motivated by the above mentioned issues 
[^, . In any case, it is gratifying that at least for solid state heterostructures 
an experimental verification of semiquantum chaos seems within reach (see 
further references in 

Our paper is organized as follows. Section II defines the model and de- 
scribes briefiy the time-dependent variational principle, from which we obtain 
our basic set of equations. In Section III we present some analytical consid- 
erations and introduce the effective semiquantum potential, which should 
help to visualize the results of our numerical studies presented in Section 
IV. In Section V we briefly compare the semiquantum dynamics with exact 
solutions of the Schrodinger equation and summarize our results. We point 
out some interesting problems related to environment-induced decoherence 
and chaotic entropy production in open quantum systems, which are left for 
future work. 

We discuss qualitatively there, how semiquantum chaos can be reconciled 
with a regular classical limit and the intrinsic linearity of quantum mechanics, 
in particular. 
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2 The Semiclassical Double- Well Oscillator 



The purpose of this section is to introduce our toy model of the simple one- 
dimensional double-well oscillator which, however, will be studied in a non- 
perturbative semiclassical approximation (TDHF, cf. below). As explained 
in the Introduction, we consider this as a necessary first step towards the 
study of a (3+l)-dimensional field theory p. The classical action of our 
model is given by 

S = j dt {\{d,^ f ~ v{v)} , (1) 

with 

v{^) ^ —f^V + , (2) 

and (p = ip{t) denotes the oscillator coordinate, which is the analogue of a 
scalar field, (p = (p{x,t), in zero space dimensions. Note that the poten- 
tial includes a negative "mass" term ~ /x^, which allows for spontaneous 
symmetry breaking in the anharmonic potential ~ (p^. Fig. 1 shows the 
potential v{ip) for the typical model parameters used throughout this work, 
/i^ = 1.0, A = 0.06 (we employ units such that h = c = 1) together with a 
typical classical phase space trajectory oscillating in one of the wells. 

The Schrodinger equation following from eq. (|l]) (coordinate representa- 
tion), 

idtij{^; t) = Hi,{^- 1) = {-1 A_ + v{v)m^- 1) , (3) 

is linear in the wave function ip. Thus, there is no mixing of its Fourier 
components, 

oo 

J dte-'-'^{if;t) , (4) 

— oo 

which allows a trivial decomposition of the full quantum dynamics w.r.t. fast 
and slow modes. This explains why at first sight one does not expect quan- 
tum chaos, no matter whether the corresponding classical system behaves 
regularly or not, and why many efforts have been made to find traces of 
classical chaos in the fully quantized system (see, e.g., Refs. fl], ^ |ll| and 
further references therein). 

We remark here that the situation for the field theory functional Schrodinger 
equation, despite its apparent similarity to eq. (^, seems even more complex. 
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since the relation of the physical (most commonly employed single-particle) 
observables to the wave functional is more complicated in this case. This is 
exemplified already by the quantum Brownian motion of a single particle in- 
teracting with a quantized electromagnetic field, phonon, or other radiation 
field i, 111, H, cf. Section V. 

Presently, we do not follow the established routes of investigations into 
quantum chaos. Rather, we study chaotic behavior generated in the semiclas- 
sical regime even for classically regular systems, such as our model defined 
by eqs. (|l|), (0). Before we derive the relevant equations of motion here, we 
pause to consider the observables of our model in more detail. Apart from 
the usual expectation values of the coordinate and its conjugate momentum, 

m = m , m = w = {-^-^) , (5) 

we will be particularly interested in expectation values of functions of ip, 

oo 

0^{t) ^ {O{0)) = I d^^{^;t)0{^)r{v;t) , (6) 



i.e. fluctuation variables. By Fourier transformation we obtain 



oo oo 



-oo — oo 



Thus, we have to find the spectrum of frequencies in 

— ip{(p;uj')ip*{^;uj' -uj) , (8) 

— oo 

which is the Fourier transform of a diagonal element of the density ma- 
trix which pertains to the pure state (coordinate representation). 
For simplicity we do not consider the off-diagonal density matrix elements 
at present, which would be needed to evaluate more general observables 
{O{0,if)). Clearly, the diagonal elements p{ip,(p]uj) yield a frequency cor- 
relation function defined at each point cp, which presents a simple and par- 
ticularly useful example out of a large variety of higher order correlation 
functions. 
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In order to illustrate the information provided by the frequency correla- 
tion function, we evaluate it formally in the full quantum case. Using the fact 
that the spectrum of the double-well Schrodinger equation corresponds 
to discrete stationary bound states, 

^n(<^;t)=0nMe''^"* , (9) 

gives 

^W(^;.;) = 27r0W(^)5(^-.;„) . (10) 
A general pure state is 

1p{ip;Uj) = 27l'^Cn(pn{'^)S{uJ - Un) ■ (H) 

n 

Then, we obtain from eqs. (|^) and (|TT|) 

p(V3, Lp;Uj) = Yj Mnn' - + (^n) , (12) 

n,n' 

i.e. a spectrum of discrete lines arising at all u = — uj'^, with a strength 

M^r^'{v) = 27rc„<,0„(<^)0:,(<^) = M:,„ . (13) 

The strength matrix M is Hermitean and time- independent, implying con- 
stant real eigenvalues. In particular, 

— oo 

p(lP, Lp;u = 0) = const ■ TrM(</)) = J dt ip{ip\t)ip*{{p;t) 

— oo 

= const ■ 27r^ |c„n0n(<^)r • (14) 

n 

Thus, the strengths of the lines depend on the actual state under consider- 
ation, i.e. the amplitudes c„, and are completely determined by the initial 
condition, ip{(f;t = 0), whereas their positions are given by the distribution 
of the energy level spacings of the system. We will employ the frequency 
correlation function as a diagnostic tool in Section IV to monitor the on- 
set of semiquantum chaos, in which case the frequency correlation function 
becomes increasingly noisy with its discrete line character eventually disap- 
pearing completely. 
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Next, we turn to the derivation of appropriate semi quantum equations 
of motion. The TDHF or Gaussian variational approximation can be easily 
obtained with the help of Dirac's variational principle The latter can be 
stated as follows: 

= , for all tfj with (M) = 1 , (15) 

dtp 

and 

' dt{mm-H)\m) , m 

i.e. requiring the effective action T defined in eq. (^6|) to be stationary 
against arbitrary variations of the normalized wave function, which vanish 
at if: ±oo, is equivalent to the exact Schrodinger equation (3). With the 
variational principle one can solve the quantum mechanical time-evolution 
problem approximately by restricting the variation of the wave function to a 
subspace of the full Hilbert space. 

In the following we work with properly normalized Gaussian trial wave 
functions, 

M^;t) ^ (27rG{t))-hxp{-[^G-\t) - za{t)][v - m? 

+^n{t)[if - , (17) 

which leads to the TDHF approximation via the evaluation of the effective 
action, eq. ([T6|) , i.e. by performing simple Gaussian integrations. Thus, we 
obtain: 

T[ijG] = Jdtl^7i^-h^-v{^) (18) 
+n[aG - 2a'G - ^-G'^ - ^v"{^)G] - ^h'v""{^)G'] , 

with v"{ip) = —fi^ + |A(^^, v""{ip) = X from eq. (0). In eq. (|T8|) we inserted 



the appropriate factors of h which exhibit the classical and semiquantum 
contributions to the effective action. Clearly, the dynamics is now described 
by the four time-dependent functions which parametrize the trial wave func- 
tion, eq. (|T^. Here cp and G play the role of the effective coordinates, while 
TT and a present the conjugate momenta. This is further illustrated by the 
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expectation values, 



{^)G = ^{t) , (-^^)G = vr(t) , (19) 



and 

{^^)g = ^'(t) + G{t) , {idt) = 7T{t)^{t) - &{t)G{t) , (20) 
which are calculated using eq. (0). Furthermore, 

((^2)G-^')^-((7rV-7r2)^ = Gi.(lG-i + 4a2G)5 

= ^(1 + (4^G')2)^>^ , (21) 

which demonstrates the uncertainty relation in the present context. Mini- 
mum uncertainty coherent states yield the lower bound on the r.h.s. of eq. 

Finally, the semiquantum equations of motion are obtained as the Eu- 
ler - Lagrange equations for the effective action, eq. (|IB]), by independent 
variations w.r.t. (f, vf, G and a: 

= 71 , (22) 

TT = -v'{^) - ^v"'{^)G , (23) 

G = AaG , (24) 

a = -2a' + ^G-' - ^v'^^) - ^v""{^)G , (25) 

with v'{ip) = — /i^<^ + |f, v"'{(f) = and v", v"" as given after eq. (p!8|). 
Several remarks are in order here: 



Equations (p2D-(p5D present a coupled set of four autonomous first- 
order nonlinear differential equations, which explains the potential for 
semiquantum chaos of the double-well oscillator, for example. 

Equivalently, by eliminating the momenta vf and a, one obtains two 
coupled second-order equations, which bear some characteristic resem- 
blance to eqs. (16) in Ref. 0, i.e. the particle production problem 
mentioned in the Introduction. 
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It follows from eq. ( pTSD that only eqs. (p2D, ( pS]) survive the simple- 
minded classical limit with h ^ 0; thus, eqs. (|2^), (p5|) present the 
dynamics of an additional effectively classical d. o. f., which arises 
through the semiclassical approximation of (Gaussian) quantum fluc- 
tuations. 



The TDHF approximation in the form of eqs. (^) - (pS]) was previously 



studied for various quantum mechanical models in Ref. ||Tj]; the potential 
for semiquantum chaos, however, has only recently been noticed based on 
a derivation via Ehrenfest's theorem 0. Note that the above equations of 
motion decouple and can be solved analytically, of course, in the harmonic 
limit, v" = const, v'" = v"" = 0. 

To conclude this section, note that the frequency correlation function, eq. 
(P), in TDHF approximation is given by 

oo 

p{ip,ip;u;)= I rfte— *(27rG)5exp{--G-Hv^-^]'} , (26) 

— CO 

using eq. (0). Clearly, the magnitude of G determines the admissible 
amount of quantum fluctuation around the classical expectation value cp, 
cf. (0), such that the integrand in eq. is not negligible. Thus, the 

relevant contributions to p{ip,ip;uj) come from those parts of a semiquan- 
tum trajectory of the system, i.e. {fit), G{t)}, which are sufficiently close to 
the point p and have a time dependence characterized approximately by the 
frequency u. 



3 Energy, Effective Potential and Instabili- 
ties 

The purpose of this section is to report some analytical considerations, which 
may help to better understand the dynamics of the semiquantum equations 



of motion, eqs. (^) - (pSj) . 

To begin with, the Hamiltonian contribution to the effective action, cf. 
eqs. (0), (|r^) together with the second of eqs. {^^, yields the conserved 
total energy for our nondissipative model. The resulting expression is still 
a rather complicated function of (p, vf, G, and a. For simplicity, we choose 
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initial conditions with zero momenta in the following, vf (t = 0) = a{t = 0) = 
0. Then, evaluating the total energy, Eq, at the minimum of the classical 
potential, for example, we obtain 

^o(G) = V^-^^ + /i^G + iAG^ • (27) 

With all other parameters kept fixed, the initial energy or, equivalently, the 
initial value G{t = 0) is the relevant control parameter for our model. 

Furthermore, comparing the initial value Eq{G) with the total energy 
evaluated at any later time, we obtain a simple measure for the numerical 
accuracy of the simulations to be discussed in Section IV. 

Secondly, in analogy to a classical system in two dimensions ("coordi- 
nates" (p and G), we consider the effective potential which is defined as the 
total minus the kinetic energy, i.e. V{lp, G) = E{ip, G; vf = a = 0), in terms 
of the total energy. Thus, 

Vi^, G) = vi^) + h[^-G-' + ^v"i^)G] + ^h\""i^)G' , (28) 

which shows the classical contribution plus quantum corrections. In Fig. 2 
we present a contour plot of equipotential hues for the effective potential. Its 
most interesting feature is the valley in the potential energy surface which 
leads to a "saddle point region" (cf. below) separating the two effective 
potential wells. Thus, our semiclassical potential shares some similarity with 
the classical "demonic potential" studied in Ref. 0. By analogy, we expect 
also in our case energy dependent transitions between regular and chaotic 
motion in the effective potential, which will be confirmed in the next section. 

Clearly, by opening the possibility to move into the semiclassical G- 
direction, our system can roll around the hill at = G = 0, i.e. around 
the equivalent of the classical local maximum of the potential at = 0. 
This "roll around" represents the quantum mechanical tunneling through 
the given potential barrier in the semiclassical TDHF approximation. Note, 
however, that the limits G — and ^ — do not commute in eq. (^) 
due to the singular term ~ /i/G, which makes the proper classical limit a 
nontrivial affair here as well as in the equations of motion, eqs. (|22|) - (^5|). 

Next, we observe that there are five fixed points of the flow described 
by eqs. (|2^) - (|25|). They correspond to the three extrema of the classical 
potential and are determined in the semiquantum case by the conditions 
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(f = TT = G = & = 0. These, by the equations of motion together with eq. 
become 7f = cr=|^ = |^ = 0, in analogy to the classical case. 



dip 

A first solution corresponding to the center of the saddle point region of 
the effective potential, cf. Fig. 2, is given by <^ = and G = /(^, j^) with 

/(a,6) =a + (6 + a3 + [62 + 26a3]5)§ + (6 + a3-[62 + 26a2]5)3 . 

For /x^ = 1.0, A = 0.06 we find G ~ 33.341. However, this turns out to be 
a small dip in the potential energy surface, since there are two fixed points 
close by. They are given by 

<^ = ±(^-3G)^ , with G = ^(1+2 cos 5) , cosa = l-(8A(-^^'^3^-^ 



A ' ' 3A ' 3' ' ' ' 3A ' ' 

For the above parameter values we find ^ ±0.106, G ~ 33.330. Finally, the 
fixed points corresponding to the minima of the two deep effective potential 
wells are given by the previous formulae for (f and G, however, replacing 
a — i> a — vr in the equation for cos a only. Thus, we obtain here (f ~ ±9.947, 
G ~ 0.355. Naturally, one may expect the trajectories to behave regularly, if 
they are initiated with sufficiently small momenta (vf , cr) sufficiently close to 
the fixed points which correspond to (local) minima of the effective potential. 

Furthermore, note that the number of fixed points has increased by two 
as compared to the classical limit. We conjecture that the number of fixed 
points would keep increasing, if we would include corrections of higher order 
in h into the effective action, eq. (plSf ), and therewith into the equations of 
motion. Equivalently, one would have to go beyond the Gaussian ansatz for 
the wave function, eq. (0), which will be further discussed in Section V. 

Finally, we turn to the investigation of instabilities arising in the equations 
of motion, which might signal the onset of chaotic behavior. It is well known 
that ordinary pertubation theory w.r.t. the nonlinear coupling, i.e. A at 
present, would not give the irregular break-up of KAM tori etc. for our 
conservative dynamical system generated by the TDHF approximation [|l]. 



11|. Therefore, we eliminate the momenta (vf , ex) from eqs. (22) - ( psf ) in 
order to rewrite them in second-order form and, thus, to expose the analogy 
to coupled (potentially unstable) oscillators. We obtain from eqs. (p2D, (|23|): 

ip + mljf^ = , (29) 



with 



mlff^^\^' + ^\G-fi' . (30) 
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From eqs. (p^, ( pSj) one finds: 



G' + ^4VG = -^p, (31) 

with 

Mf^^ = + AG - 2/^2 . (32) 

Obviously, both oscillators can become unstable due to exponentially growing 
modes for a finite range of (initial) values for and a. Both effective mass 
squared terms become negative for 

< G < ^ - . (33) 
A 

Note that the repulsive term on the r.h.s. of eq. (|3lD keeps G > at all 
times. Thus, there is a region in the positive (^-G-halfplane, where we expect 
trajectories with slightly different initial conditions (keeping for simplicity 
7f = (^ = cr = G = initially) to diverge exponentially. 

However, eq. (|33|) seems to provide at best a criterion for the local sensi- 
tivity to initial conditions. In particular, for initial conditions corresponding 
to high energy, i.e. sufficiently small G{t = 0) according to eq. (^), one 
may expect the trajectories to behave rather regularly, since the unstable 
region in the (^-G-plane by eq. (p3D becomes a smaller and smaller fraction 
of the accessible one with increasing energy. These effects will be seen in our 
numerical results to be presented in the following section. 



4 Numerical Results 

In this section we present the results of numerically studying the solutions 



of the semiquantum equations of motion, eqs. ( Ua) - (|25D , for a variety of 
initial conditions. However, we always fixed vf = a = initially and chose 
G = 0.1 together with various values for to represent initial configurations 
with varying total energy, E{7t = a = 0) = V{i^,G), cf. eq. (^). Generally, 
we did not explore the full set of all initial conditions at a given energy, which 
would require much more extensive numerical calculations. 

To begin with, we show in Fig. 3 the results for two energies which are 
very close to the absolute minimum, Emin ~ —24.2948; Emin is obtained at 
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the minima of the effective potential ((^ = ±9.947, G = 0.355) discussed 
previously. Thus, we observe in the left column of Fig. 3, which pertains 
to the lowest energy, a rather simple regular motion characterized by the 
basic oscillator frequency ujq = \pl and a second frequency at twice this 
value. Note that the value of c^o is obtained by expanding the classical 
potential, eq. (|^), around one of its minima and identifying the coefficient 
of the harmonic term, i.e. hJ^ = /i^ = 1.0, at present. Together with the 



frequency correlation function, cf. eqs. (|]), p^), and ([Ty|), we present two 
Poincare sections in the phase space planes of the two conjugate pairs of 
variables, {vr, (^} and {a, G}, respectively. 

Furthermore, a full trajectory is portrayed in the last row in Fig. 3. 
Clearly, a slight increase of the energy, with the corresponding results pre- 
sented in the second column here, immediately leads to the admixture of 
higher multiples of the basic frequency cuq- In this case, the Poincare sec- 
tions still look fairly simple, but the picture of a trajectory already shows 
a rather convoluted structure. Obviously, the increasing complexity of tra- 
jectories can be expected to be reflected in the more and more detailed fine 
structure of the frequency correlation function as the energy is increased 
further. 

Before we proceed, a few remarks are in order here concerning the numer- 
ical solution of the equations of motion and the evaluation of the frequency 
correlation function, in particular. 

We have used Mathematica ll5| to solve eqs. (|^) - (^) numerically. In 
particular, we employed the ND solve function which iteratively and adap- 
tively solves the equations to arbitrary accuracy (up to the working precision 
of the computer, of course) which we chose to be 10 digits in the final solu- 
tion. The solutions for various initial conditions were calculated for varying 
lengths of time depending on the complexity of the trajectory. Thus, the 
results in Fig. 3 were obtained with 200 time units {h = c = 1) while one of 
the Poincare sections in Fig. 4 (see below) was obtained from a solution of 
length 10,000 time units (the longest of our trajectories). 

As mentioned earlier, one measure of the accuracy of our solution is to 
compare the total energy E at the beginning and ending of each trajectory. 
For E = —1.2661 initially and the initial conditions stated above, we obtain 
E = -1.2660 after 10,000 time units. 

Evaluating the numerical solutions to eqs. ( p2D - (|^) at discrete time in- 
tervals, we find the frequency correlation function by discrete Fourier trans- 
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form, i.e. by replacing the integral in eq. ( p6D with a sum over the discrete 
time intervals of the trajectory. For our purposes it was sufficient to take 
the discrete time interval At = 0.5 which is more or less arbitrary. It is 
sufficient in the sense that the discrete nature of p{ip, ip; u) is clearly evident, 
or not, depending on the particular value of the total energy and the initial 
conditions. Tests with smaller values of At on selected trajectories produced 
the same low frequency spectrum. 

The length of the trajectories, r, used to calculate p{ip, ip; u) was also 
somewhat arbitrary, r must be long enough to capture the basic low fre- 
quency structure. Again, our objective was to observe the qualitative behav- 
ior of the frequency correlation function rather than make a precise determi- 
nation of its content, and we found that trajectories of length 200 time units 
were sufficient for this purpose. For the regular trajectories, increasing r 
did not alter the frequency spectrum, and for the trajectories that exhibited 
chaotic character, increasing the trajectory length did not change the overall 
character of the frequency correlation function. 

Finally, we note that p{p, p; u) as calculated above is normalized to 
times its continuum value, eq. (^), where r = NAt. 

To construct Poincare sections, we simply evaluated each solution for a 
particular total energy at an arbitrary fixed value of one of the four variables 
while the others were allowed to vary. In practice, a small variation in the 
fixed d.o.f. was allowed to reduce the number of evaluations of the solution, 
i.e. our two-dimensional slice of phase space had a small but finite thick- 
ness. Also, to observe the characteristic phase space filling of the chaotic 
trajectories, r was increased to as much as 10,000 time units. 

Next, we present in Fig. 4 characteristic Poincare sections for various 
higher energies. At E = —12.44 (initialized with G = 0.01) we still observe 
regular behavior, which is also obvious from the frequency correlation func- 
tion (not shown). Notice that the trajectory is still completely confined to 
the right effective potential well, cf. Fig. 2, where it was initialized. Clearly, 
one expects that different initial conditions at the same energy would shift 
this Poincare section w.r.t. the coordinate axes and distort it as well. At 
E = —6.61 the structure of the trajectories becomes much more complicated 
with the Poincare section showing indications of break-up of KAM tori. This 
is better represented in the frequency correlation function, see. Fig. 5, where 
there are still outstanding discrete lines. However, they appear already con- 
siderably shifted w.r.t. to (integer multiples of) uq- Furthermore, there is 
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now a broad low-level background at fractional frequencies corresponding to 
the rather spiky character of the trajectory itself. From the Poincare section 
it is obvious that the trajectory begins to run into the valley of the effective 
potential, cf. Fig. 2, which is the semiclassical equivalent to (attempts at) 
quantum mechanical tunneling to the other potential well < 0). 

At E = —1.266 and E = 15.6 the Poincare sections in Fig. 4 both show 
the "stochastic sea outlining islands of stability" which is the typical feature 
of irregular motion resulting from the break-up of KAM tori |[l|, |Tl| . Note that 
the lower energy here still corresponds to trajectories which classically could 
not get to the other side of the barrier at = (height f (0) = 0), i.e. they 
have to tunnel there. At E = 15.6 trajectories run well above the classical 
barrier and would yield regular motion in the classical limit. However, both 
cases show semiquantum chaos as a result of quantum corrections. 

Finally, in Fig. 5 the frequency correlation function is also shown for two 
other energies, E = —3.686 and E = 201.8, which reveal further interesting 
features of our system. First, as expected, in the high energy limit the behav- 
ior becomes regular again, since the effective potential essentially becomes a 
single well with a perturbation of decreasing influence at the center. This is 
obvious from the appearance of the frequency correlation function. Presum- 
ably, its spectrum could be related to the quantum mechanical energy levels 
of the anharmonic (^"^-oscillator relevant here at high excitation energy. 

Secondly, and quite surprisingly indeed, we find at E = —3.686 regular 
behavior with a large number of discrete lines in the spectrum, which are 
rather evenly spaced and of comparable strength. It seems to originate from 
a resonance phenomenon leading to a fairly simple closed trajectory confined 
to one well, see Fig. 6. In distinction to the cases of Fig. 3, where the motion 
results from the quasi linear superposition of a few harmonic frequencies, at 
present the trajectory appears simple despite being fully nonlinear. This 
orderly behavior arises at an energy which lies in between cases with much 
more irregular if not fully chaotic trajectories, cf. Fig. 5. We did not 
yet investigate in full detail the sensitivity of this phenomenon to changes 
of the energy control parameter etc. In particular, it is an open question, 
whether other regular energy values like this one exist within the irregular 
intermediate energy regime. Preliminary investigations, however, show that 
this "resonating trajectory" is even sensitive to where precisely it is initiated 
on an equipotential contour, cf. Fig. 2. Presumably, we encounter a regime 
with stable KAM tori side by side with irregular domains |]1|, |IT|, which is 
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already indicated by the Poincare sections in Fig. 4. 

The results presented here naturally raise a number of more general ques- 
tions. They will be discussed in the following section. 

5 Conclusions 

Presently, we have been studying semiquantum chaos in a simple one-dimensional 
double-well oscillator model. In the Introduction we explained the motivation 
for this and related studies by field theory considerations and nonrelativistic 
many-body problems 0, |, ||, |^. In this context our present work presents 
first steps towards a more detailed understanding of those more complicated 
systems and their behavior in the semiclassical regime, in particular. 

Our work is based on the simple observation that a semiclassical descrip- 
tion (employing the Gaussian variational principle or the TDHF approxi- 
mation in our case, cf. Section III) necessarily involves the introduction of 
additional degrees of freedom beyond the classical ones, which may define the 
model. The corresponding additional coordinate and conjugate momentum 
variables are related to quantum fluctuations of the classical ones, i.e. they 
represent correlation functions of varying order (cf. below) in the quantum 
mechanical sense 0. The classical and nonclassical degrees of freedom are 
governed by a coupled set of first-order autonomous differential equations, 
cf. eqs. (^) - (p5l). They are necessarily highly nonlinear, which to a major 
extent is the case independently of the particular dynamical model. Herein 
we find the potential for deterministic semiquantum chaos in this effective 
nondissipative Hamiltonian system. 

Our numerical results presented in Section IV confirm the above consid- 
erations, in general. In Section III we argued that the total energy in our 
model, as determined, for example, by the parameters of the initial state 
wave function, cf. eq. (|T7p, serves as the control parameter. Depending on 
the energy we find regular and chaotic trajectories as solutions of the semi- 
classical equations of motion. Besides Poincare sections we mainly employed 
a suitably defined frequency correlation function, cf. eqs. (H), (p6D , which is 
related to the density matrix of the system to monitor its behavior. Here the 
transition to chaos could be expected to show up as a drowning of the discrete 
line structure deduced in Section II in a rising broadband noise. This was 
clearly observed in the numerical results, see, e.g.. Fig. 5. However, quite 
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surprisingly we found that the behavior of the system as a function of the 
energy control parameter changes in a rather complex way. Apparently, reg- 
ular and irregular regimes alternate with each other, where only the simplest 
cases (low/high energy limit) can be easily understood. The simple-minded 
stability analysis in Section III allows to outline qualitatively some regime of 
irregular behavior, but is quantitatively not sufficient to point out more of 
the detailed structure as a function of energy. 

Clearly, our numerical studies show that there is room for more surprises 
in the behavior of the system depending on the control parameter. Fur- 
thermore, we did not fully explore the phase space of our system at a given 
energy: Are all initial conditions equivalent in that the corresponding tra- 
jectories are either regular or irregular? These problems, being intimately 
related to the question of ergodicity in the semiclassical model, can only 
be addressed through more extensive numerical work, which is beyond the 
scope of our present paper. Similarly, we did not investigate the nature of 
the transition to chaos in this system, nor did we look into details of the 
control parameter dependence of the frequency correlation function spectra. 

However, our results together with the above considerations also point 
towards several issues of more general interest, which we discuss in the re- 
maining part of this section. 

It is convenient to introduce a formal representation of the exact wave 
function and Schrodinger equation for any quantum mechanical system which 
generalizes the Gaussian ansatz and variational equations of motion, eqs. 
(p^ and p^ ) - (p5D, repectively. We rewrite the wave function, 

oo 

^p{x, t) = exp /(x, t) = exp ^ . (34) 

n=0 

This ansatz, of course, can only be considered as formal, since we do not 
know how to work out the normalization condition, 

1 = J dxtp*tp = J dx exp{2Re/(x,t)} = F[Re/„(t)] , (35) 

in general. Nevertheless, pretending we know the highly nonlinear functional 
F, we obtain a representation of the exact Schrodinger equation, cf. eq. (|^), 
as an infinite autonomous set of coupled first-order differential equations, 

^/o = -^/2 + ^^ , (36) 
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(39) 



(38) 



(37) 



(n + 2)(n + l) 
2 



fn+2 1 



(41) 



where we assumed for definiteness a simple polynomial potential, v{x) = 
ax + bx"^ + cx^ + dx^. Several observations seem worth mentioning here: 

• The above exact representation of the full quantum dynamics is not 
limited to one dimension and can be extended to field theory using 
the functional Schrodinger picture, cf. |^, for example. However, 
the power series ansatz, eq. (|3^ , is not general enough or the most 
convenient one to cover all interesting cases and is used here only for 
heuristic reasons. 

• There is an essential nonlinearity due to the normalization condition 
which, however, appears only in the first of the infinite set of coupled 
equations. 

• Different polynomial type interactions enter the equations as different 
constant parameters only. 

Based on these observations we conjecture the following. Any finite trun- 
cation of sufficiently high order [n < N) in such a nonlinear expansion scheme 
for the full quantum dynamics, even of a classically regular system, may lead 
to semiquantum chaos (in the sense presented in this paper) depending on the 
energy control parameter. In general, we expect additional classical conserva- 
tion laws to lead to a higher dimensional control parameter space. Since the 
one and only non-linearity originates from the wave function normalization 
condition, cf. eqs. (pSj), (|36|), and enters the flow equations in a completely 
different manner than any interactions, we expect some still to be specified 
universality of the phenomena of semiquantum chaos. 

Finally, knowing that any exact wave function solution of the linear 
Schrodinger equation cannot itself behave chaotically as discussed in Section 
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II, we are led to speculate that the infinite number of linearly coupled cor- 
relation degrees of freedom represented by eqs. (|36D - (^Tj) effectively restore 
the linear behavior of the full quantum dynamics. The latter here was only 
undermined by a single nonlinear constraint entering as a scale independent 
nonlinear coupling into one of the flow equations! 

This brings us to the flnal topic of our discussion. In view of the linear 
character of the Schrodinger wave equation and particularly for classically 
regular systems one might be tempted to consider semiquantum chaos as 
a spurious effect caused by a truncation of the full quantum dynamics as 
described above or any other bad approximation thereof. 

First of all, to illustrate the quality of the time- dependent Hartree-Fock 
approximation we include here Fig. 7 showing some typical comparisons be- 
tween TDHF results and those of exact solutions of the Schrodinger equation 
for given initial conditions. Overall the expectation value of the position of 
the wave packet, for example, is always well described for sufficiently short 
times. Here the scale should be set by the characteristic time for the quan- 
tum mechanical spreading of the initial wave packet. Clearly, this has to 
depend on the energy at present. Even the long-time behavior is qualitively 
reproduced rather well in Fig. 7, no matter whether the TDHF trajectories 
are regular or chaotic at a particular energy. There is one notable excep- 
tion here at E = —1.266, where TDHF for the chosen initial conditions is 
fully tunneling by a "roll-around" as discussed in Section II, whereas the 
major part of the exact wave packet stays is the initial well for long times. 
Obviously, TDHF misses some of the interference structure of full quantum 
mechanics. For further comparisons of TDHF and exact results in different 



quantum mechanical models we refer to Ref. [14 



In any case, however, one may argue that for a closed system the full 
quantum description is always better than any approximation. Thus, the 
way out of this apparent no-go situation, in order to establish the physical 
relevance of semiquantum chaos, is to allow for dissipative effects. Necessar- 
ily, one has to consider the interaction of the quantum system under study 
with its environment . In this way, the observation of semiquantum chaos is 
tied to the identification of suitable open quantum systems which through 
interactions with their decohering environment [16| are driven into the semi- 
classical regime. 

Various steps in approaching a physically relevant dynamically enforced 



classical limit (not naively setting h ^ 0) have been made P, O, |17|, ITS 
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ranging in context from cosmology through strong interactions, questions of 
foundations of statistical mechanics, and the interpretation of quantum me- 
chanics. Presently, with the examples discussed in the Introduction in mind, 
it seems most relevant to find systems, e.g. in solid state physics ("quan- 
tum dots" [l^), where environment- induced quantum decoherence and con- 



sequently semiquantum chaos become experimentally accessible. 

Based on the example treated in Ref. @], a quantum particle showing 
essentially classical behavior due to its interaction with a suitable quantized 
radiation field, the following model suggests itself for further study: the 
double-well oscillator governing a quantum particle coupled to a decohering 
radiation field. Before any detailed results, it would be interesting to see 
the effect of a deliberately chosen environment on the relevant set of flow 
equations, such as eqs. (|36|) - (^1]) above. We plan to address some of these 
fascinating problems in a future publication. 
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Figure Captions: 



Fig. 1: The double-well oscillator potential, eq. with fi"^ = 1.0, A = 0.06. 
Also shown is a typical classical phase space trajectory (0 = (^). 

Fig. 2: The effective semiquantum potential, eq. (|28|) , for = <^ > (V is 
symmetric in 0). The inset shows the detailed shape of the contours 
close to the positive minimum (i.e. fixed point) of V. 

Fig. 3: For two low energies (left/right column) the following TDHF results 
are shown (top to bottom): frequency correlation function, Poincare 
sections in the two different coordinate/conjugate momentum planes, 
full trajectory (see text for details). 

Fig. 4: Poincare sections at various values of the energy (control parameter) 
illustrating the transition from regular to chaotic motion. 

Fig. 5: Frequency correlation functions at various energies. Semiquantum chaos 
shows up as broadband noise in the expected line spectrum. There is 
an unexpected resonance effect at E = —3.68601. 

Fig. 6: The nonlinear resonance trajectory discussed in the text. 

Fig. 7: Comparisons between exact solutions of the time-dependent Schrodinger 
equation (full lines) and TDHF results (symbols) at various energies. 
Shown here is the expectation value of the position of the particle, 
which is initially a Gaussian wave packet. 
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